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Abstract 



It is commonly accepted that the peak effect (PE) in the critical current density of type II 
superconductors is a consequence of an order-disorder transition in the vortex lattice (VL). Ex- 
amination of vortex lattice configurations (VLCs) in its vicinity requires the use of experimental 
techniques that exclude current induced VL reorganization. By means of linear ac susceptibility 
experiments in the Campbell regime, where vortices are forced to oscillate (harmonically) around 
their effective pinning potentials, we explore quasi-static stable and metastable VLCs in NbSe2 
single crystals near the PE. We identify three different regions: for T < Ti(H), stable VLCs are 
maximally ordered. For T > T 2 (H) configurations are fully disordered and no metastability is 
observed. In the T\ < T < Ti region we find temperature dependent stable configurations with 
intermediate degree of disorder, possibly associated to coexistence of ordered and disordered 
lattices throughout the PE. A simple estimation of the equilibrium proportion of ordered and 
disordered domains is provided. 



The behavior of systems under the influence of both, thermal and quenched disorder has been 
an issue of intense research over the last 50 years in different areas of condensed matter physics. 
Quenched disorder can broaden first order phase transitions inducing phenomena as phase separa- 
tion (under active study in the context of manganites for instance) and if sufficiently strong can 
even change completely the properties of this system producing the appearance of new phases (as 
it is the case in many glassy states) [L. 
*email:pasquini@df.uba.ar 
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In this context, the vortex lattice (VL) in superconductors provides models systems where elastic 
and pinning interactions together with thermal fluctuations compete. A fingerprints of complex 
behavior in this system is the Peak Effect (PE), an anomalous non- monotonous dependence of the 
critical current density J c with both temperature and magnetic field [2] . The origin and nature of 
this behavior are still controversial issues. 

Phenomenological pictures based in an order-disorder (O-D) transition from a quasi-ordered 
Bragg glass (BG) [3 j to a disordered phase with proliferation of topological defects [I] explain a 
broad amount of related experimental results. However, the underlying physics of the transition, 
the structural symmetry of the disordered phase |5] and its connection with a clean limit melting 
transition [7J |6] are issues that remain unsolved. While increasing amount of experimental evi- 
dence [8] suggests that the nature of the PE could depend on the various materials, the marked 
metastability and history effects are common facts reported both in low [9j [101 E] and in high T c 
|12] materials. 

In traditional superconductors, neutron diffraction experiments show a clear change in the struc- 
ture factor lattice at the PE indicating an O-D transition [T3]. In NbSe2, metastable coexisting 
regions with different J c have been directly observed [T3]; this fact as well as transport measure- 
ments in the Corbino geometry [TT] and experiments of magnetization assisted by a shaking ac field 
|10] , suggest a first order transition where the ac magnetic field assists an equilibration process from 
a supercooled disordered metastable phase to an ordered stable BG phase. 

PE and history effects are often studied by means of transport experiments. In standard 
transport techniques, vortex undergo some type of current-induced reorganization and therefore 
the original vortex lattice configuration (VLC) is not accessible. To overcome this problem, Xiao et 
al. [7] have developed a sophisticated ultrafast transport technique and have shown the existence 
of an enlarged crossing phase boundary between ordered and disordered phases in NbSe2 crystals. 
A well known and appropriate technique to measure the pinned VL response without external 
disturbances is ac susceptibility restricted to the linear Campbell regime [15] where a very small 
ac field h a superimposed to the dc field H is applied, forcing vortices to perform small (harmonic) 
oscillations inside their effective pinning potential wells. 

In this letter we explore quasi-static VLCs in the Campbell regime in NbSe2 single crystals 
in the vicinity of the PE following different thermal, magnetic and dynamical histories. Essential 
to our work, we avoid any measurement induced VL reorganization. We are able to access the 
corresponding stable configuration at each temperature (field) by applying a shaking ac field. For 
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the first time, the pinning potential curvature of stable quasi-static VLCs is measured throughout 
the PE, resulting in a clear identification of low temperature maximally ordered, high temperature 
maximally disordered and stable configurations with intermediate degree of disorder at interme- 
diate temperatures. In the last case, the stable VLCs are separated by energy barriers, possibly 
related with surfaces between coexisting ordered and disordered VL domains. We can also ob- 
serve the evolution of the pinning potential with temperature in metastable VLCs and confirm the 
spontaneous ordering in cooling processes |16j . 

Results shown here correspond to a NbSe2 single crystal [T7] of approximate dimensions (0.5 x 
0.5 x 0.03) mm 3 , with Tc = 7.30 K (defined as the mid point of the ac susceptibility linear 
transition at H = 0) and ATc = ±0.02 K. Ac susceptibility has been measured with a home made 
susceptometer based in the mutual inductance technique. All fields are parallel to the c axis of the 
sample. 

In the Campbell limit, the inductive component of ac susceptibility x' 1S determined by the 
geometry and the dimensionless parameter Ar/D (Ar is the real penetration depth and D is a 
characteristic sample dimension). In this regime the imaginary penetration depth Aj « Ar . The 
curvature of the effective pinning potential well is the Labusch constant ai that can be numerically 
estimated from Ar = {A 2 l +4>qB / 'Altai) 1 / 2 ; where A^ is the London penetration depth (for references 
and details of the numerical procedure see Ref . [T8] ) . 

Fig. 1 shows typical X"(T) (panel a) and X'(T') (panel b) field cool cooling (FCC) and field cool 
warming (FCW) curves in the linear (h a = 0.025 Oe, black) and non-linear regime {h a = 0.32 Oe, 
gray) at H = 320 Oe. The temperatures at which x' shows a maximum or a minimum (T° n and 
T p in Ref. [19]) slightly depend on h a or history . Above a temperature T2 > T p FCW and FCC 
curves merge. Below Ti , A/ < 0.1 Ar so dissipation is very small (x" ~ 0) an d X' 1S frequency 
independent in the kHz range, indicating a Campbell regime. In the non-linear regime, where 
vortices are forced out of their pinning sites, both FCC and FCW curves display PE whereas in 
the linear Campbell regime it appears only in the warming process. As the PE is a signature of 
an O-D transition, its absence in the FCC Campbell regime indicates that the VL nucleates at 
high T and remains trapped in a metastable disordered and strongly pinned configuration. On 
the other hand, even if ac induced currents are not able to assist the VL, the system warms with 
lower pinning potential (larger penetration depth). In this framework, the hysteresis displayed in 
the Campbell regime must imply some kind of spontaneous re-ordering at low temperatures. This 
ordering during FCC was recently observed by using a time resolved transport technique |16j . In 
the inset of Fig. 2, X'(T) recorded in various warming processes from different low temperatures is 
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Figure 1: Typical X"(T) (panel a) and X'(T) (panel b) FCC and FCW curves in the linear (black curves) 
and non-linear regime (gray curves) measured at /= 30 kHz. Arrows indicate the direction of temper- 
ature variation. T\ and T2 limit the T regions of stable ordered and disordered phases (see text). Inset: 
Linear FCW curves at various dc fields ploted as a function of (T — Tc(H)). An abrupt PE develops for 
H > 200 Oe. The zero dc field transition is also shown. 

shown. A lower initial temperature results in a less pinned (more ordered) VLC. We have tested 
that the warming curves are independent on the cooling/warming rate, and are identical if the 
cooling process is performed without measurement. 

This spontaneous ordering at low temperatures, can be enhanced by applying a large shaking 
ac field. This shaking field allows the system to explore the free-energy landscape and to reach a 
more ordered and stable VLC. This feature is illustrated in the main panel of Fig. 2. The result of 
the FCW process from 4.2 K shown in the inset (gray line curve in both panels ) is compared with 
various warming curves (symbols) recorded after shaking the lattice at different temperatures T s h 
(large vertical arrows). The following process has been performed: the system has been measured 
in the linear regime during FCC procedure and stabilized at T^. Then the measurement was 
interrupted and a sinusoidal shaking ac field of 3.2 Oe - 30 kHz was applied during 30 s. The 
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Figure 2: (Color online). Inset: X'(T) recorded in various warming processes from different T : there is 
a progressive spontaneous ordering in the cooling process. Main panel: FCC and FCW (4.2 K) processes 
(gray lines) compared with curves (various symbols) recorded before and after shaking the VL at different 
temperatures T s h (large vertical arrows). Vertical (magenta) line remarks the huge x' change after shaking. 
All shaken VLCs collapse to the same reversible curve. Small arrows indicate the change in temperature 
with time. 

measurement {h a = 25 mOe) was resumed in a warming procedure starting from T s h leading to 
a drastic decrease of pinning (vertical color lines in the figure). The various symbols are hard to 
distinguish because all warming procedures collapse to the same curve. This curve has remarkable 
characteristics below a temperature T\ < T° n (dotted line in the figure): first, it is reversible; 
second, after any additional shaking, with any waveform or amplitude, the system remains in this 
configuration [20]. These features constitute a first important result: below T\ the stable VLCs 
(attained via smooth temperature changes) are continuously connected in a minimum of the free- 
energy landscape. These are the less pinned (more ordered) VLCs at each T. This picture is 
consistent with a stable ordered BG below T\ where changes in T only produce elastic (reversible) 
deformations. 

The description changes completely when T% < T < T^: the more ordered VLCs become 
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unstable and the PE develops. Again, the system does not reach spontaneously a stable VLC: 
warming (cooling) leads to metastable over ordered (disordered) configurations. In every case, a 
large shaking field assists the system in accessing a stable VLC. The qualitative difference observed 
above T\ is that, once a stable VLC has been reached, any smooth temperature variation results 
in new metastable states: the stable configuration at the new temperature is only reached through 
shaking; consequently, no reversibility is possible. This remarkable feature is one of the central 
issues of this Letter, and is illustrated in Fig. 3 where x' values in the linear regime corresponding 
to the stable VLCs at various T are plotted in large symbols. The complete experiment shown 
in Fig. 3 is the following: First a warming process (small black dots) was performed stabilizing 
temperature successively at very small intervals AT^ [2T]. At each T s h (indicated by black large 
vertical arrows) the measurement was interrupted and the same shaking protocol described above 
was performed, in way to reach the stable VLC (open triangles). At fixed temperature, additional 
perturbations do not modify the response and the system remains in this stable configuration. Then 
the measurement on warming was resumed (small black dots) until the next T s h was reached. As 
T increases the corresponding stable VLC is more disordered. 

Once T > Ti was reached, a cooling process (small gray dots) was performed and a similar 
protocol was followed, shaking the system at the same temperatures (gray vertical arrows). Notice 
that x'{T s h) obtained after shaking in warming (hollow triangles) and cooling (gray circles) is the 
same (within experimental resolution). This central result supports the stability of VLCs with 
intermediate order, and rules out the possibility that they arise from a harder metastability or slow 
dynamics. The more ordered (warming curve after shaking at T s ^ < T\) and the more disordered 
FCC curve are shown for reference in black line. 

To quantify the pinning strength on the various VLCs, the effective Labush parameter «l(T) 
corresponding to each VLC has been estimated from the x' values in the linear regime at 90 kHz in 
the region where Xj « Xr. The result is shown in Fig. 4a. Black curves correspond to metastable 
VLCs evolving in FCC and FCW processes whereas the gray curve corresponds to ordered VLCs 
attained by shaking at 6.4 K. Hollow circles show ai corresponding to stable VLCs between T\ 
and T2 extracted from Fig. 3. It can be seen that all the ai(T) curves approach at low T. At 
T > T\ the ordered VLCs (gray curve) become metastable. The ai corresponding to the stable 
VLC increases more than five times in the PE, indicating a huge increase of disorder, qualitatively 
different to the smooth and monotonically increasing ai(T) observed in YBCO crystals |8| 20J. The 
ordered (reversible) gray curve below T± corresponds to a BG, free of dislocations, and the black 
(reversible) curve above T2 corresponds to a completely disordered phase. The stable character 
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Figure 3: Warming (W) and cooling (C) experiment between T± and T2. Vertical down (up) arrows 
identify the T s h in the W (C) process. At each T s h the lattice is shaken to find the stable VLC (large open 
triangles in W and large gray circles in C ). The stable VLC warms (cools) in a metastable VLC (small 
dots). The more ordered and disordered curves are also shown in black line for reference. 

of the VLCs with intermediate degree of disorder for T\ < T < T2, together with additional 
evidence of coexisting domains in NbSe2 p3] , allows us to propose a scenario of equilibrium phase 
separation. In this framework, the metastability would arise from an exceeding proportion of 
ordered or disordered domains. To give a rough estimate of the proportion of ordered phase, we 
proceeded in the following way: we extrapolated the ordered Labusch parameter a or d(T) to zero 
at Tc(H) (dotted line in Fig. 4a) and we proposed otdi s (T) = 7 a or d(T) to fit the black curve above 
T2 (dash-dotted line in Fig. 4a). We then used the calculated ai{T) to estimate the proportion of 
ordered phase n or d(T) = (otdi s (T) — cxl{T)) / (otdi s {T) — a or d(T)) corresponding to each VLC. The 
result is shown in Fig. 4b. 

We arrive then to the following scenario. There are three different temperature regions along 
the PE: for T < T\(H), stable VLCs are maximally ordered and accessible by the application of 
a shaking ac field. We identify these ordered stable VLCs as a BG, continuously connected by 
elastic deformations in a minimum of the free-energy landscape. The robust metaestable states 
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Figure 4: (a) a£,(T) in log. scale corresponding to the various stable and metastable VLCs. otord {ddis) is 
a fitting curve that extrapolates a(T) in a completely ordered (disordered) VL. (b) Estimated proportion 
of ordered phase in each VLC. Arrows indicate the change in temperature with time 

below T\ show spontaneous ordering at low T as elastic interactions overcome pinning forces, al- 
though remain separated from the BG by a finite barrier in all the experiemntal range, down to 
4.2K. For T > T2{H) configurations are fully disordered and no metastability is observed. We 
identify T2{H) with Ts(H), the spinodal line in Ref [TH].The evolution from the ordered BG to 
the high temperature disordered phase, in the small T\ < T < T2 region, occurs by a gradual 
increase in the proportion of disorder. The equilibrium VLCs are only accessed from over-ordered 
or over-disordered configurations by a shaking ac field, that provides sufficient energy for plastic 
(irreversible) deformations to access the stable VLC. Plastic changes imply the creation or annihila- 
tion of dislocations, probably related to the movement of domain walls in a scenario of equilibrium 
phase separation. In conclusion, by measuring for the first time the pinning potential curvature of 
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stable quasi-static VLCs, avoiding vortex lattice reorganization, we are able to present a consistent 
picture that describes the vortex lattice physics in the PE region in NbSe2 single crystals. 
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